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Abstract 

The surface detector array of the Pierre Auger Observatory consists of 
1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) 
generated by ultra-high-energy cosmic rays. We describe the trigger hierar- 
chy, from the identification of candidate showers at the level of a single detec- 
tor, amongst a large background (mainly random single cosmic ray muons), 
up to the selection of real events and the rejection of random coincidences. 
Such trigger makes the surface detector array fully efficient for the detection 
of EAS with energy above 3 x 10^^ eV, for all zenith angles between 0° and 
60°, independently of the position of the impact point and of the mass of 
the primary particle. In these range of energies and angles, the exposure of 
the surface array can be determined purely on the basis of the geometrical 
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acceptance. 



Key words: Ultra high energy cosmic rays, Auger Observatory, Extensive 
air showers. Trigger, Exposure 
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1 1. Introduction 

2 The main objective of the Pierre Auger Collaboration is to measure the 

3 flux, arrival direction distribution and mass composition of cosmic rays from 

4 fti 10^^ eV up to the highest energies. Due to the very low fluxes at these 

5 energies, cosmic rays have to be measured through the extensive air showers 

6 (EAS) they produce in the atmosphere. 

7 The Pierre Auger Observatory, located near Malargiie, Argentina, at 1400 

8 m asl, detects EAS in two independent and complementary ways. It includes 

9 a surface detector array (SD), consisting of 1600 water-Cherenkov detectors 

10 [l| on a triangular grid of 1.5 km spacing covering an area of approximately 

11 3000 km^, which detects the secondary particles at ground level and thus 

12 samples their lateral density distribution. The surface detector array is over- 
is looked by a fluorescence detector (FD) consisting of 24 telescopes at four 

14 sites, which measure the fluorescence light emitted along the path of the air- 

15 showers and thus traces their longitudinal development [2l| . Showers detected 

16 by both detectors are called hybrid events and they are characterised more 

17 accurately with respect to direction and energy than using either technique 

18 alone. However, the livetime of the FD is limited to ^ 13%, as it only op- 

19 erates on clear, moonless nights [2]. The bulk of data is provided by the SD 

20 with its nearly 100% livetime. The study of the trigger and the determina- 

21 tion of the aperture of the SD is thus essential for the physics aims of the 

22 Pierre Auger Observatory. 

23 The SD data acquisition (DAQ) trigger must fulfill both physical and 

24 technical requirements. The main limitation to the rate of recordable events 

25 comes from the wireless communication system which connects the surface 

26 detectors to the central campus. The latter must serve continuously 1600 

27 stations spread over 3000 km^, each using an emitter consuming less than 

28 1 W power to transmit to collectors as far as 40 km away. The maximum 

29 sustainable rate of events per detector is less than one per hour, to be com- 

30 pared to the 3 kHz counting rate per station, due to the atmospheric muon 

31 flux. The trigger thus must reduce the single station rate, without induc- 
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32 ing loss of physics events. It must also allow data acquisition down to the 

33 lowest possible energy. To deal with all these requirements, the design of 

34 the DAQ SD trigger (described in section 3) has been realised in a hierar- 

35 chical form, where at each level the single station rate becomes less and less, 

36 by means of discrimination against background stricter and stricter. At the 

37 same time, the DAQ trigger is designed to allow the storage of the largest 

38 possible number of EAS candidates . 

39 The ultimate discrimination of EAS from chance events due to combina- 

40 torial coincidences among the surface detectors is performed off-line through 

41 a selection of physics events, and of detectors participating in each of them. 

42 The event selection procedure is hierarchical too, it is described in section 4. 

43 In section 5.1, we show that the trigger and event selection hierarchy 

44 makes the array fully efficient for the detection of showers above 3 x 10^^ eV. 

45 We restrict ourselves to this energy range for the calculation of the exposure 

46 (described in section 5.2), which is simply proportional to the observation 

47 time and to the geometrical size of the SD array. Under these conditions the 

48 calculation of the exposure is very robust and almost devoid of systematic 

49 uncertainties. Therefore it is straightforward to calculate the cosmic ray flux 

50 as the ratio of the number of collected events to the effective, as it was done 

51 in the measurement of the cosmic ray spectrum by the surface detector of 

52 Auger Q. 



53 2. The surface detector of the Pierre Auger Observatory 

54 Each water Cherenkov detector of the surface array has a 10 m^ water 

55 surface area and 1.2 m water depth, with three 9" photomultiplier tubes 

56 (PMTs) looking through optical coupling material into the water volume, 

57 which is contained in a Tyvek® reflective liner |l|, Each detector operates 

58 autonomously, with its own electronics and communications systems powered 

59 by solar energy. Each PMT provides two signals, which are digitised by 40 

60 MHz 10-bit Flash Analog to Digital Converters (FADCs). One signal is 

61 directly taken from the anode of the PMT, and the other signal is provided 

62 by the last dynode, amplifled and inverted within the PMT base electronics to 

63 a total signal nominally 32 times the anode signal. The two signals are used 

64 to provide sufficient dynamic range to cover with good precision both the 

65 signals produced in the detectors near the shower core (~ 1000 particles/^s) 

66 and those produced far from the shower core (~ 1 particle/yus). Each FADC 

67 bin corresponds to 25 ns 
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68 The signals from the three PMTs are sent to a central data acquisition system 

69 (CDAS) once a candidate shower event triggers the surface detector array (see 

70 section 3.2). The total bandwidth available for data transmission from the 

71 detectors to the CDAS is 1200 bits per second, which precludes the possibility 

72 of any remote calibration. For this reason, the calibration of each detector 

73 is performed locally and automatically. It relies on the measurement of the 

74 average charge collected by a PMT from the Cherenkov light produced by a 

75 vertical and central through-going muon, Qvem 0. The water-Cherenkov 

76 detector in its normal configuration has no way to select only vertical muons. 

77 However, the distribution of the light of atmospheric muons produces a peak 

78 in the charge distribution, Qy^'^M (or VEM in short), as well as a peak in 

79 that of the pulse height, Iy'^\^, both of them being proportional to those 

80 produced by a vertical through-going muon. The calibration parameters are 

81 determined with 2% accuracy every 60 s and returned to the CDAS with 

82 each event. Due to the limited bandwidth, the first level triggers are also 

83 performed locally. These triggers (section 3.1) are set in electronic units 

84 (channels): the reference unit is ly'^M- 

85 With respect to shower reconstruction, the signals recorded by the detec- 

86 tors - evaluated by integrating the FADC bins of the traces - are converted 

87 to units of Qvem- These are fitted with a measured Lateral Distribution 



88 Function (LDF) [ll|, that describes S'(r), the signals as a function of dis- 

89 tance r from the shower core, to find the signal at 1000 m, S'(IOOO) (6|. The 

90 variation of S'(IOOO) with zenith angle 9 arising from the evolution of the 

91 shower, is quantified by applying the constant integral intensity cut method 

92 0, justified by the approximately isotropic flux of primary cosmic rays. An 

93 energy estimator for each event, independent of 9, is 5*38, the S'(IOOO) that 

94 EAS would have produced had they arrived at the median zenith angle, 38°. 

95 The energy corresponding to each 5*38 is then obtained through a calibration 

96 with the fluorescence detector based on a subset of high-quality hybrid events 
1. 



97 



98 3. The DAQ trigger system of the surface detector array 

99 The trigger for the surface detector array is hierarchical. Two levels of 

100 trigger (called Tl and T2) are formed at each detector. T2 triggers are com- 

101 bined with those from other detectors and examined for spatial and temporal 

102 correlations, leading to an array trigger (T3). The T3 trigger initiates data 
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103 acquisition and storage. The logic of this trigger system is summarised in 

104 figure dJ 



Tl (Single Station Level) 



T2 (Single Station Level) 



T3 (CD AS level) 
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1 bin 
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Figure 1: Schematics of the hierarchy of the trigger system of the Auger surface detector. 
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3.1. Single detector triggers 

The Tl triggers data acquisition in each water Cherenkov detector: data 
are stored on the local disk for 10 s waiting for a possible T3. Two indepen- 
dent trigger modes are implemented as Tl, having been conceived to detect, 
in a complementary way, the electromagnetic and muonic components of an 
air-shower. The first Tl mode is a simple threshold trigger (TH) which re- 
quires the coincidence of the three PMTs each above 1.75 /{f^yll- This trigger 
is used to select large signals that are not necessarily spread in time. It is par- 
ticularly effective for the detection of very inclined showers that have crossed 
a large amount of atmosphere and are consequently dominantly muonic. The 
TH-Tl trigger is used to reduce the rate due to atmospheric muons from Ri3 
kHz to ^100 Hz. The second Tl mode makes use of the fact that, for other 
than very inclined showers or signals from more vertical showers very close 
to the shower axis, the arrival of particles and photons at the detector is 
dispersed in time [g, 0|. For example, at 1000 m from the axis of a vertical 
shower, the time for the signal from a water-Cherenkov detector to rise from 
10 to 50% is about 300 ns. The second mode is designated the "Time-over- 
Threshold" trigger (ToT) and at least 13 bins (i.e. >325 ns) in 120 FADC 



^For detectors with only two (one) operating PMTs the threshold is 2 (2.8) 1^'^°-'^ 



'VEM- 
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123 bins of a sliding window of 3/is are required to be above a threshold of 0.2 

124 Iy'eIi coincidence in 2 out of 3 PMT^. This trigger is intended to select 

125 sequences of small signals spread in time. The ToT trigger is thus optimised 

126 for the detection of near-by, low energy showers, dominated by the electro- 

127 magnetic component, or for high-energy showers where the core is distant. 

128 The time spread arises from a combination of scattering (electromagnetic 

129 component) and geometrical effects (muons) as discussed in js], |9| where de- 

130 tails are given of how the time spread depends on distance and zenith angle. 

131 Since the average signal duration of a single muon is only about 150 ns, the 

132 time spread of the ToT (325 ns) is very efficient at eliminating the random 

133 muon background. The ToT rate at each detector is < 2Hz and is mainly 

134 due to the occurrence of two muons arriving within 3/is, the duration of the 

135 sliding window. 

136 The T2 is applied in the station controller to reduce to about 20 Hz the 

137 rate of events per detector. This reduction is done to cope with the band- 

138 width of the communication system between the detectors and the central 

139 campus. The T2 triggers, namely their time stamp and the kind of T2, are 

140 sent to the CDAS for the formation of the trigger of the array. All ToT-Tl 

141 triggers are promoted to the T2 level, whereas TH-Tl triggers are requested 

142 to pass a further higher threshold of 3.2 Iy'em coincidence among the three 

143 PMTtl. The rates of the TH-T2 triggers are rather uniform in the detectors 

144 over the whole array within a few percent, while those due to the ToT-T2 

145 are less uniform. This is due to the fact that the ToT is very sensitive to the 

146 shape of the signal, this in turn depending on the characteristics of the water, 

147 the reflective liner in the detector and the electronic pulse shaper. However, 

148 the lack of uniformity of the trigger response over the array does not af- 

149 feet the event selection or reconstruction above the energy corresponding to 

150 saturated acceptance. 

151 3.2. Trigger of the surface array 

152 The third level trigger, T3, initiates the central data acquisition from the 

153 array. It is formed at the CDAS, and it is based on the spatial and temporal 

154 combination of T2. Once a T3 is formed, all FADC signals from detectors 



^For detectors with only two (one) operating PMTs, the algorithm is applied to two 
(one) PMTs. 

•^For detectors with only two (one) operating PMTs the threshold is set to 3.8 (4.5) 

jpeak 
^VEM- 
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155 passing the T2 are sent to the CDAS, as well as those from detectors passing 

156 the Tl but not the T2, provided that they are within 30 /is of the T3. 

157 The trigger of the array is realised in two modes. The first T3 mode 

158 requires the coincidence of at least three detectors that have passed the ToT 

159 condition and that meet the requirement of a minimum of compactness, 

160 namely, one of the detectors must have one of its closest neighbours and 

161 one of its second closest neighbours triggered. It is called "T0T2C1&3C2", 

162 where C„ indicates the n*'' set of neighbours (see figure [2]). Once the spatial 

163 coincidence is verified, timing criteria are imposed: each T2 must be within 

164 (6 + 5C„)/xs of the first one. An example of such T3 configuration is shown 

165 in figure [21 left. Since the ToT as a local trigger has very low background, 

166 this trigger selects predominantly physics events. The rate of this T3 with 

167 the full array in operation is around 1600 events per day, meaning that each 

168 detector participates in an event about 3 times per day. This trigger is 

169 extremely pure since 90% of the selected events are real showers and it is 

170 mostly efficient for showers below 60°. The 10% remaining are caused by 

171 chance coincidences due to the permissive timing criteria. The second T3 

172 mode is more permissive. It requires a four-fold coincidence of any T2 with 

173 a moderate compactness. Namely, among the four fired detectors, within 

174 appropriate time windows, at least one must be in the first set of neighbours 

175 from a selected station (Ci), another one must be in the second set (C2) 

176 and the last one can be as far as in the fourth set (C4). This trigger is 

177 called "2C1&3C2&4C4". Concerning timing criteria, we apply the same logic 

178 as for the "T0T2C1&3C2". An example of such T3 configuration, is shown 

179 in figure [2l right. Such a trigger is efficient for the detection of horizontal 

180 showers that, being rich in muons, generate in the detectors signals that have 

181 a narrow time spread, with triggered detectors having wide-spread patterns 

182 on the ground. With the full array configuration, this trigger selects about 

183 1200 events per day, out of which about 10% are real showers. 
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Figure 2: Example of T3 configurations: the 3-fold T3 mode ToT2Ci&3C2 is shown on 
the left and the 4-fold mode 2Ci&;3C2&4C4 on the right (see text for the definitions). CI, 
C2, C3, C4 indicate the first, second, third and fourth sets of neighbours, respectively at 
1.5, 3, 4.5 and 6 km from a given detector. 

184 3.3. Efficiency of the single detector trigger 

185 The single detector trigger probability as a function of the signal, V{S), 

186 besides being important for the determination of the efficiency of the trigger 

187 of the array, is also of use in the event reconstruction where non-triggered 

188 detectors are included up to 10 km from a triggered one |10|. 

189 The Tl efficiency versus signal in the detector, 'P{S), is determined by 

190 using the very large statistics of EAS (~ 10^) recorded by the surface detector 

191 array. For each detected EAS, and each participating detector, we measure 

192 the trigger probability V{S) as the ratio in different bins of 9 and 

193 5(1000), of the number of triggered stations, N^, to the total number of 

194 active stations, Nqn- S is the expected signal at a detector, based upon the 

195 LDF fitted from the measured values from each detector, and 5(1000) is the 

196 signal strength at 1 km, as derived from this fit. Since V{S) is obtained from 

197 events that actually produced a T3, the method is biased by events with a 

198 positive fiuctuation in the signal. This bias can be corrected by Monte Carlo 

199 simulations and is found to be negligible at energies above around 3 x 10^^ 

200 eV. Limiting the analysis to showers with 5*38 > 16 VEM (corresponding to 

201 about 3 X lO^^eV), the trigger probability versus signal is derived averaging 

202 over all the bins in 6 and S(IOOO). This is shown in figure |3] (circles): the 

203 probability becomes > 0.95% for 5 ~ 10 VEM. This result is confirmed by an 

204 independent analysis that makes use of showers triggering certain detectors 

205 that have been specially located very close to one another. The surface 

206 array has seven positions in which three detectors (so called triplets) have 

207 been deployed at 11 m from each other. In each triplet, only one detector 
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208 (master) sends T2 to CD AS, while the other two (slaves) are independently 

209 read out each time a T3 is generated and if they pass the Tl. For each 

210 slave, the trigger probability versus recorded signal S is derived from the 

211 ratio between the number of events where both slaves have triggered and 

212 the number of events where only the other one has triggered. Depending if 

213 one or two slaves have triggered, S is either the signal of the only triggered 

214 detector or the average of the two. 
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Figure 3: Single detector trigger probability as a function of the signal in the detector, 
ViS), obtained from triplets data (triangles) and from showers data with E > 3 x 10^* eV 
(circles) . 

215 From the analysis of about 10000 events, and combining the probabilities 

216 for the two slaves, V{S) is obtained and it is shown in figure [2] (triangles), in 

217 good agreement with the one obtained by showers data. 



4. Event selection of the surface detector array for showers with 
zenith angle below 60° 

A selection of physics events and of detectors belonging to each event is 
made after data acquisition. Indeed, a large number of chance coincidence 
events is expected due to the large number of possible combinations among 
the single detectors. We focus here on the selection of events between 0° and 
60° since more inclined showers have different properties and require specific 
selection criteria described elsewhere 12 . 
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226 Two successive levels of selection are implemented. The first one (physics 

227 trigger) is based on space and time configurations of the detector, besides 

228 taking into account the kind of trigger in each of them. The second one 

229 (fiducial trigger) requires that the shower selected by the physics trigger is 

230 contained within the array boundaries, to guarantee the accuracy of the event 

231 reconstruction both in terms of arrival direction and energy determination. 

232 The logic of this off-line trigger system and its connection to the DAQ triggers 

233 is summarised in figure |H 




Figure 4: Schematics of the hierarchy of the event selection of the Auger surface detector. 

234 Physics trigger 

235 The physics trigger, T4, is needed to select real showers from the set of 

236 stored T3 data. Two criteria are defined, with different aims. The first T4 

237 criterion, so-called 3ToT, requires 3 nearby stations, passing the T2-ToT, in 

238 a triangular pattern. It requires additionally that the times of the signals in 

239 the 3 stations fit to a plane shower front moving at the speed of the light. 

240 The number of chance coincidence passing the 3 ToT condition over the full 

241 array is less than one per day, thanks to the very low rate of the T2-ToT. 

242 Due to their compactness, events with zenith angles below 60°are selected 

243 with high efficiency, i. e. more than 98%. 

244 The second T4 criterion, so called 4C1, requires 4 nearby stations, with 

245 no condition on the kind of T2. In this case also, it is required that the 

246 times of the signals in the 4 stations fit to a plane shower front moving at 
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247 the speed of the hght. This 4C1 trigger brings to ^100% the efficiency for 

248 showers below 60°. 

249 The zenith angle distribution of events selected by the T4 criteria is shown 

250 in figure O left, in the unfilled histogram for 3ToT, and in the filled one for 

251 the 4C1 that are not 3ToT: the two criteria are clearly complementary, the 

252 latter favouring the selection of events with larger zenith angles. In figure El 

253 right, the energy distributions of events selected by the two different criteria 

254 are shown: those selected by 3ToT have a median energy around 6 x 10^^ 

255 eV, while for those selected by 4C1 it is around 3 x 10^® eV. 
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Figure 5: Left: Angular (left) and energy (right) distribution of events selected by the T4 
triggers: 3ToT (unfilled histogram), and 4C1, not ToT (filled histogram). 

256 Besides disentangling accidental events, there is also the need to identify, 

257 and reject, accidental detectors in real events, i.e. detectors whose signals 

258 are by chance in time with the others, but that in fact are not part of the 

259 event. To this aim, we define a "seed" made by 3 neighbouring detectors in 

260 a non-aligned configuration. If there is more than one triangle of stations, 

261 the seed with the highest total signal is chosen. If the T4 is a 3 ToT, only 

262 ToT detectors can be considered to define the seed; if it is a 4C1, also TH 

263 detectors can be included. Once the triangle has been determined, the arrival 

264 direction is estimated by fitting the arrival times of the signals to a plane 

265 shower front moving with the speed of light. Subsequently, all other detectors 

266 are examined, and are defined as accidental if their time delay with respect 

267 to the front plane is outside a time window of [—2 yus : +1 /is]. Detectors that 
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268 have no triggered neighbours within 3 km are always removed. 

269 After the selection chain (both event selection and accidental detectors 

270 removal), 99.9% of the selected events pass the full reconstruction procedure, 

271 that is arrival direction, core position and 5(1000) are determined. 

272 4-^- Fiducial trigger 

273 The need for a fiducial trigger, T5, mainly arises from events falling close 

274 to the border of the array, where a part of the shower may be missing. In 

275 figure [6] a hybrid event is shown, that triggered the SD and one of the FD 

276 telescopes, where a part of the SD information is missing due to its position 

277 on the border of the array. 




Figure 6: Example of a hybrid, non-T5, event: the event falls on the border of the SD 
array, triggering only four detectors. Filled circles indicate the triggered ones, open circles 
the non-triggered active ones. The dimensions of the filled circles are proportional to 
the measured signal. The shower detector plane reconstructed by FD (dash-dotted line) 
indicates that the core is within the triangle of detectors. The SD only reconstruction 
places it outside the array (cross), artificially increasing the event energy. 

278 Such events could have wrong core positions, and consequently, incor- 

279 rect energies, as in this example where the energy derived by SD is more 

280 than 4 times larger than the one estimated by FD (1.4 x 10^^ eV instead 

281 of 3 X 10^^ eV). The main task of the fiducial trigger is thus to select only 

282 events well contained in the array, ensuring that the shower core is properly 

283 reconstructed. 
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The fiducial trigger should be applied a priori on the events, to be inde- 
pendent of the reconstruction procedure. The T5 adopted requires that the 
detector with the highest signal has all its 6 closest neighbours working at 
the time of the event (i.e., it must be surrounded by a working hexagon). 
This ensures adequate containment of the event inside the array. Even in the 
case of a high energy event that falls inside, but close to the border of the 
array, where part of the data may be missing, information from the seven 
detectors closest to the shower core ensures a proper reconstruction. Apply- 
ing this condition, the maximum statistical uncertainty in the reconstructed 



S(IOOO) due to event sampling by the array is ^ 3% |10|. It has to be noted 
that this criterion also discards events that, though contained, fall close to 
a non-working detector: this is an important issue because, due to the large 
number of detectors distributed over 3000 km^, about 1% of the detectors are 
expected to be not functioning at any moment, even with constant detector 
maintenance. For the fully completed array, and taking this into account, 
the application of the T5 condition reduces the effective area by 10% with 
respect to the nominal one. 

Finally, the use of the fiducial trigger allows the effective area of the array 
to saturate to the geometrical one above a certain primary energy. Indeed, 
with no conditions on event containment, the acceptance would increase with 
increasing energy, since showers falling outside the borders of the array might 
still trigger sufficient detectors to be recorded; the higher their energy, the 
farther the distance. 



307 5. Aperture and exposure of the surface detector array for showers 

308 vi^ith zenith angle below 60 degrees 

309 The aperture of the surface detector array is given by the effective area 

310 integrated over solid angle. When the trigger and event selection have full 

311 efficiency, i.e. when the acceptance does not depend on the nature of the 

312 primary particle, its energy or arrival direction, the effective area coincides 

313 with the geometrical one. In subsection 5.1, the energy above which the 

314 acceptance saturates is derived. In section 5.2, the calculation of the exposure 

315 above this energy is detailed. 

316 5.1. Determination of the acceptance saturation energy 

317 I. Prom SD data. The acceptance saturation energy, Esat, is deter- 

318 mined using two different methods that use events recorded by the surface 
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detector array. In the first one, starting from detected showers, mock events 
are generated by fluctuating the amplitude of the signals recorded in each 



detector and their arrival time. Such fluctuations are measured 13|, [iJl by 
using twin detectors located at 11 m from each other. To each simulated 
event, the full trigger and event selection chain are applied. From the ratio 
of the number of triggered events to the simulated, the trigger efficiency is 
obtained as a function of energy, as shown in figure [7] (triangles) . As can 
be seen, the trigger probability becomes almost unity (> 97%) at energy 
~ 3 X 10^^ eV for all angles between 0° and 60°. The fact that the method 
is based on the use of showers that actually triggered the array may bias 
the estimation of the trigger probability at low energy. However, it does not 
bias the result on the trigger probability close to full efficiency, and hence on 

EsAT- 
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Figure 7: Trigger efficiency as a function of energy, derived from SD data (triangles) and 
hybrid data (circles). 



332 II. From hybrid data. The hybrid data sample is composed of events 

333 observed by the FD and that triggered at least one SD detector: consequently, 

334 it has an intrinsically lower energy threshold than the SD. For each bin in 

335 energy (of width 0.2 in log^o(E)), the number of events that pass the SD 
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336 trigger out of the total number of events are counted. To avoid biases from 



337 primary composition, the same data selection criteria as in [15| are used. 

338 Additionally, in analogy with the T5, to avoid the effects of the borders 

339 of the array, it is required that the detector used in the hybrid geometry 

340 reconstruction is surrounded by 6 active detectors. The trigger efficiency of 

341 the surface detector array is found to be saturated (> 97% ) for energies 

342 above 3 x 10^^ eV, as shown in figure [7] (circles), in agreement with what is 

343 obtained by the analysis of SD data alone. 

344 III. Cross-check with simulations. Esat is finally cross-checked us- 

345 ing full shower and detector simulations. The simulation sample consists of 

346 about 5000 proton, 5000 photon and 3000 iron showers simulated using COR- 



347 SIKA [16| with zenith angle distributed as sin 9 cos 9 (9 <60°) and energies 

348 ranging between 10^'' eV and 10^^'^ eV in steps of 0.25 (0.5 for photons) in 

349 logio(E). The showers are generated using QGSJET-II and FLUKA [isj 

350 for high and low energy hadronic interactions, respectively. Core positions 

351 are uniformly distributed at ground and each shower is used five times, each 

352 time with a different core position, to increase the statistics with a negligible 

353 degree of correlation. The surface detector array response is simulated using 

354 Geant4 [19] within the framework provided by the Off line software [20]. The 

355 resulting trigger probability as a function of the Monte Carlo energy for pro- 

356 ton, iron and photon primaries is shown in Figure [8] for 0° < < 60°. Due 

357 to their larger muon content, at low energies iron primaries are slightly more 

358 efficient at triggering the array than protons. However, the trigger becomes 

359 fully efficient at 3 x 10^*^ eV, both for proton and iron primaries, in different 

360 intervals of zenith angles. It is important to notice that the trigger efficiency 

361 for photons is much lower. This is because photons tend to produce deeper 

362 showers that are poor in muons. 
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Figure 8: SD trigger efficiency as a function of Monte Carlo energy E for proton (circles), 
iron (triangles) and photon primaries (squares) and zenith angle integrated up to 60°. 
Lines are drawn only to guide the eyes. 



363 5.2. Calculation of the integrated exposure 

364 The studies described above have shown that the full efficiency of the 

365 SD trigger and event selection is reached at 3 x 10^^ eV. Above this energy, 

366 the calculation of the exposure is based solely on the determination of the 

367 geometrical aperture and of the observation time. 

368 With respect to the aperture, the choice of a fiducial trigger based on 

369 hexagons, as explained in section 14. 2^ allows us to exploit the regularity of 

370 the array very simply. The aperture of the array is obtained as a multiple of 

371 the aperture of an elemental hexagon cell, Oceii, defined as any active detector 

372 with six active neighbours, as shown in figure [9l 
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Figure 9: Scheme of an hexagon of detectors: the elemental hexagon cell, flccii, is the 
shaded area around the central detector. 

373 At full efficiency, the detection area per cell is 1.95 km^. The correspond- 

374 ing aperture for showers with 9 < 60° is then Oceii — 4.59 km^ sr. The number 

375 of cells, Nceii{t), is not constant over time due to temporary problems at the 

376 detectors (e. g. failures of electronics, power supply, communication system, 

377 etc.). Nccii{t) is monitored second by second: we show in figure ITU] the evo- 

378 lution of Nccii{t) between the start of the data taking, January 2004, and 

379 December 2008. Such precise monitoring of the array configurations allows 

380 us to exploit data during all deployment phases, clearly visible in the figure, 

381 as well as during unstable periods as during, for example, January 2008 when 

382 huge storms affected the communication system. 
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Figure 10: Evolution of the number of hexagonal cells (see text) between January 1**, 
2004 and December 31"*, 2008 



383 The second-by-second monitoring provides at the same time the aperture 

384 of the array per second, Oceii x A^ceii(^); as well as the observation time with 

385 high precision. To calculate the integrated exposure over a given period of 

386 time, the aperture of the array, Nceii{t) x aceii, is integrated over the number of 

387 live seconds. This calculation is expected to be very precise, since it is based 

388 on a purely geometrical aperture and a very good time precision. However 

389 both the determination of A^ceii(i^) and of the observation time are affected 

390 by uncertainties. 

391 Concerning the determination of Nccii{t), to evaluate the uncertainty in 

392 the number of active detectors, a check of the consistency of the event rate of 

393 each detector with its running time, determined from the monitoring system, 

394 is performed. The uncertainty derived from this study is added to that due 

395 to errors of communication between the station and the DAQ, which are also 

396 monitored. Overall, the uncertainty on the determination of Ncc\\{t) amounts 

397 to about 1.5%. 

398 For the determination of the observation time, and related uncertainty, 

399 the dead time that is unaccounted for in the second by second monitoring 
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400 of the array, is taken into accounto. To determine these, an empirical tech- 

401 nique is exploited, based on the study of the distribution of the arrival times 

402 of events, under the reasonable hypothesis that they follow a Poisson dis- 

403 tribution. Given the constant rate A for the T5 event rate per hexagon, 

404 A ~ 1.4 X 10^^ event per second per hexagon, the probability P that the 

405 time interval T between two consecutive T5 events be larger than T is given 

406 by: P{T) = e~^'^ . We define intervals as dead time if the Poisson probability 

407 of their occurrance is less than 10~^. As an example, we show in figure [11] 

408 the distribution of time differences for events acquired in 2008. The distribu- 

409 tion is exponential with a time constant of 72.4 seconds, as expected for the 

410 above value of A and the observed average number of live hexagons during 

411 that year. In the figure, the points outside the filled area show those time 

412 intervals that have occurred with a Poisson probability less than 10~^. 
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Figure 11: Distribution of time differences between events in 2008. The points outside the 
filled area show the dead times (see text). The exponential fit is shown as a dashed line 
in the inset where the histograms are zoomed. 



413 The identified dead times generally correspond to periods of software 



'''This dead time can be due either to problems in the communication between the 
stations and the CDAS or to problems of data storage in the stations 
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414 modifications at the level either of the single detectors or of the CDAS. These 

415 were rather frequent during the deployment phase of the surface detector 

416 array, which lasted until June 2008. The uncertainty in the determination 

417 of the livetime is estimated to be around 1%. Between January 2004 and 

418 December 2008, the livetime of the surface detector array data acquisition is 

419 96%. Hidden dead times reduce the effective livetime to 87%, the reduction 

420 being mostly due to the two first years of operation. However, due to the 

421 growth of the surface detector array, their impact on the total integrated 

422 exposure is a reduction of only 3%. 

423 6. Conclusions 

424 The DAQ trigger of the surface detector array of the Pierre Auger Ob- 

425 servatory is organised in a hierarchical way, starting at the level of the single 

426 detector (Tl, T2) up to the data acquisition (T3). The selection of events 

427 below 60° takes place off-line, and it is also hierarchical (T4, T5). The whole 

428 chain, from the single detector trigger, up to event selection, is able to re- 

429 duce the counting rate of the single detector from about 3 kHz, due mainly 

430 to single, uncorrelated, cosmic muons, down to about 3 x 10~^ Hz. This 

431 final rate is due to extensive air showers, more than 99% of which pass the 

432 reconstruction chain. 

433 In spite of the large number of detectors and the possible number of 

434 chance events due to combinatorial coincidences among the detectors, the 

435 high-purity Time Over Threshold trigger enables the main trigger of the 

436 array to be kept at the level of a 3-fold coincidence, thus extending the 

437 range of physics that can be studied. Such a trigger, together with the 

438 event selection strategy, allows the acquisition and reconstruction of about 

439 one cosmic ray shower per minute, with median energy around 6 x 10^^ eV. 

440 Moreover, it makes the surface detector array fully efficient for showers due to 

441 primary cosmic rays above 3 x 10^® eV, independent of their mass and arrival 

442 directions. The trigger provides at the same time a larger overlapping energy 

443 region with the FD, which is naturally efficient at lower ener gies , allowing 

444 the measurement of the cosmic ray spectrum down to 10^* eV |2l|. 




445 Above 3 X 10^^ eV, the calculation of the exposure is purely geometrical, 

446 being the integration of the geometrical aperture over the observation time. 

447 Both of them are known with high precision, so that the overall uncertainty 

448 on the integrated exposure is less than 3%. The integrated SD exposure as a 

449 function of time is shown in figure [T2| from January 2004 to December 2008: 
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450 at the end of the period it amounts to 12790 ±380 km^ sr yr. Even though the 

451 SD was under continuous deployment until June 2008, the effective livetime of 

452 the surface detector array averaged over all the five years is high, being 87%. 

453 The effective livetime of the SD is 96% for 2008 alone: with this livetime and 

454 the full surface detector array deployed, the exposure is expected to increase 

455 by about 500 km^ sr yr per month. 




Date 

Figure 12: Evolution of the integrated exposure between January l'**, 2004 and December 
3r\ 2008 
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